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We consider an electronic phase separation process that generates regions of different charge 
densities, or local dopings, as the origin of the inhomogeneous charge density of high T c supercon- 
ductors. We show that it gives rise to a phase boundary potential between such doping disordered 
regions or grains. The Bogliubov-deGennes self-consistent calculations in this disordered medium 
yield position dependent superconducting gaps which are, for all dopings, smaller than those de- 
rived from the local density of states with a pseudogap behavior. Studying these two sets of gaps for 
different temperatures and dopings, we are able to reproduce many many observed properties of su- 
perconducting cuprates. This scenario is consistent with a resistivity transition driven by Josephson 
coupling among the superconducting grains. 



PACS numbers: 74.20.-z, 74.25.Dw, 74.62.En, 74.72. Kf 

The origin of the superconducting gap associated with 
the superconducting state and the relation to the pseu- 
dogap above the transition temperature (T c ) remains one 
of the central questions in high-T c research. In this let- 
ter, we show that the solution of this important problem 
is connected to an electronic granular structure derived 
from an electronic phase separation (EPS) transition. 

There are steadily accumulating evidences that the 
charge distribution is microscopically inhomogeneous in 
the C11O2 planes of high temperature superconductors 
(HTSC)[l|-Q. Low temperature Scanning Tunneling Mi- 
croscopy (STM) has demonstrated non-uniform energy 
gaps A that vary on the length scale of nmQ over 
the whole surface of these materials. These gaps have 
two types of shapes [3], and some remain well above 
the superconducting critical temperature T c (p)[|], 
More recently, the existence of two types of energy 
gaps have been observed on electronic Raman scattering 
experiments @, STM data^-Q, Angle Resolved Photon 
Emission (ARPES)[I3 and combined STM-ARPES[n|. 
However, the origin and even the existence of this two- 
gap picture is still a matter of debate [13, [l^. 

A possible reason to this complex behavior is an 
EPS transition driven by the minimization of the free 
energy [14}. This minimization is due to the formation 
of anti-ferromagnetic (AF) regions of almost zero local 
doping [pi ss 0) with lower free energies than the re- 
gions with large local doping level (j>i « 2p), where p 
is the average hole doping of the sample [14|. As p in- 
creases, the Coulomb repulsion in the large local doping 
regions generates a high energy cost to the phase separa- 
tion process and the EPS ceases in the overdoped region, 
in agreement with the disappearance of the local AF fluc- 
tuations [l5j] . This cooperative phenomenon gives rise 
to (grain) boundary potentials in the C11O2 planes, sin- 
gle particle bound states and intra-grain superconductiv- 
ity at low temperature. These superconducting regions 
develop Josephson couplings among them and the resis- 



tivity transition T c (p) occurs when the Josephson energy 
Ej(p) is equal to K B T c (p). 

We used the theory of Gahn-HiUiard (CH)[H-[l| to 
binary alloys to describe the EPS transition. The transi- 
tion order parameter is the difference between the local 
and the average charge density u(p, i, T) = (p(i, T)—p)/p. 
In Fig. ([T]), we show a typical density map with the two 
(hole-rich and hole-poor) main solutions given by differ- 
ent colors. 




FIG. 1: (color online) The density map simulation of the 
inhomogeneous charge density on a 100 x 100 sites. In the 
inset, the GL potential with the two minima representing the 
two (high and low density) solutions [3] and the potential 
barrier V g t between them. Marked white points are located 
where some calculations shown here were made. 



In the inset of Fig. ([T]) we plot the usual potential from 
the Ginzburg-Landau (GL) free energy expansion used 
in the CH equation [14| with its two minima and the po- 
tential barrier V g b = V g b(p,T) between them. The barri- 
ers generate shallow potential wells in the CuOi planes 
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that always have an exponential small bound statep^j. 
These bound states lower the kinetic energy to allow, un- 
der favorable conditions, local hole pair formation in the 
antiadiabatic limit [20|. 

The calculated u(i,T) (or p(i,T)) density map, as 
shown in Fig.©, is used as the initial input and 
it is maintained fixed throughout the self-consistent 
Bogoliubov-deGennes (BdG) calculations. Notice that 
the charge inhomogeneity p(i,T and V g b(p, T) are cor- 
related. In what follows, the parameters involved in 
V g b(p, T) are chosen to match the average local density of 
states (LDOS) gaps measured by low temperature STM 
on 0.11 < p < 0.19 Bi2212 compounds [|]. All others 
parameters are similar to values previously used[2~ll. |22T|. 

The effect of the temperature in the potential is taken 
into account using V g b(T) ~ (1 — (T/Tpg) 15 , as demon- 
strated by Cahn and Hilliard[lj|. Thus we were able 
to obtain the intra-grain or local superconducting tem- 
peratures T c (i), i.e., the onset temperature for the d- 
wave superconducting gap Ad(i, T) at a given location i. 
The largest value of all T c (i) determines the temperature 
T on (p) that marks the onset of local superconductivity of 
the sample. 

Similar to granular superconductors |23| . the super- 
conducting transition occurs in two steps: first by the 
appearing of intra-grain superconductivity and than by 
Josephson coupling with phase locking at a lower temper- 
ature. This approach provides a clear interpretation to 
the superconducting amplitude and the measured quasi- 
particles dispersion above Tr- (p) fl2l. Il3l| . 

By using the theory of granular superconductors [HI to 
these electronic grains and the calculated average super- 
conducting amplitudes Af (T,p) = £)f A d (T, i,p)/N, 
where N is the total number of sites, we can estimate the 
values of T c (p). 
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Where A™(T, p) is the average of the local supercon- 
ducting gaps A d (i,p,T) on a N x N (N = 28, 36 and 
42) square lattice which are plotted in the top panel of 
Fig-©. The inset shows the schematically single par- 
ticle levels at a shallow puddle whose walls are propor- 
tional to V g b- R n is the normal resistance of a given com- 
pound, which is proportional to the planar resistivity p a b 
measurements [25( on the La2~ v Sr p CuOi series. In the 
low panel of Fig.©, the Josephson coupling Ej(p,T) is 
plotted together with the thermal energy ksT whose in- 
tersection yields the critical temperature T c (p), as shown 
in the inset. The values are in reasonable agreement with 
the Bi2212 T c (p), as expected, since V g b was chosen to 
match the Bi2212 low temperature LDOS[3[. 

In the BdG approach, the symmetric local density of 
states (LDOS) is proportional to the spectral function[26] 
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FIG. 2: Top panel, the average Ad(p), and in the inset, the 
schematic single particle and superconducting energy levels 
in a grain. In the low panel, the thermal energy ksT and the 
Josephson coupling among superconducting grains Ej(p,T) 
for some selected doping values as function of T. The inter- 
sections give T c (p), as plotted in the inset. 



and may be written as 

Ni(T,V gb ,eV) = ^[|u„(x l )| 2 + |«„(x J )| 2 ] x 

[f' n (eV-E n ) + f' n (eV + E n )]. (2) 

The prime is the derivative with respect to the argu- 
ment. u n ,v n and E n are respectively the eigenvectors 
and eigenvalues of the BdG matrix equationpj], l2lj . 
/„ is the Fermi function and V is the applied volt- 
age. Ni(T, Vgb, eV) = LDOS{Vgb) is proportional to 
the tunneling conductance dl/dV, and we probe the ef- 
fects of the inhomogeneous field by examining the ra- 
tio LDOS{V gb ^ 0)/LDOS(V gb = 0). LDOS{V gb = 
0) contains the inhomogeneous charge distribution and 
LDOS(V g b 7^ 0) vanishes around the Fermi energy due 
to the superconducting gap and the single particle lev- 
els in the grains. This LDOS ratio yields well-defined 
peaks and converges to the unity at large bias, in close 
agreement to similar LDOS ratios calculated from STM 
measurements @ . 

We calculated the LDOS at many values of p and T 
to compare with the STM 0-0] and ARPESflMl data > 
but here we present detailed results only for an under- 
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doped (p — 0.11) and an overdoped (p = 0.21) sample. 
Thus, in Fig. ^ we show the results at two representa- 
tive places (shown as white dots in the middle of Fig. (1)) 
of the p = 0.11 sample at various temperatures. In this 
case, the low temperature LDOS display very few well 
defined (coherent) peaks but the gaps are similar to the 
STM data of McElroy et a! 3 . The low temperature 
superconducting gap A^ produces small anomalies that 
are marked by arrows, as reported by some STM data@- 
This behavior was already noticed by our previous 
workflij]. Then the gaps derived from the LDOS peaks 
are identified with the local pseudogap Apo(i,p,T) be- 
cause they are larger than the local Ad(i,p,T) and be- 
cause they vanish at a higher temperature T*(i) than 
To(i). 

Let us investigate in detail the relationship between 
Apc{i,p, T) and Ad(i,p, T). For the p — 0.11 compound 
with T c (0.11) « 65K at a typical hole-poor place with 
Pi = 0.02, we obtain a low temperature Apg ~ 80meV 
and Ad = 12meV. In a representative hole-rich loca- 
tion with p(i) w 0.23 there is an interchange of in- 
tensity: Ape ~ 50meV and the superconducting gap 
Ad{T = 0) = 33meV. We observe that this inversion 
in the values of Apq and A^ is common to all com- 
pounds. For this p = 0.11 sample, at all sites, Apg(j) 
and Ad(i) vanish with increasing temperature almost to- 
gether near T = 145K. At the pi = 0.02 hole poor grain, 
A PG and A d vanish near T*(i) = T c (i) = 147K, and at 
the pi = 0.21 hole rich grain, they close at T*(i) = 148 
and T c (i) — 146K. We see that, for underdoped sam- 
ples, T*(i,p) T c (i,p) and they remain much above the 
resistivity transition T c (0.11) w 65K, i.e., both gaps are 
present in the pseudogap phase up to T* . 

The difference between T*(i,p) and T c (i,p) increases 
continuously with increasing doping. Calculations with 
p = 0.16 show that in some grains T*(i) ps T c (i) + 15K 
while in some other locations T*(i) » T c (i). By in- 
creasing the doping p, this difference increases as shown 
in Fig. (0| for p = 0.21. For this compound with 
T c (0.21) « 72K, most of the LDOS present well de- 
fined coherent peaks as observed by McElroy et al||. 
On the top panel, we show the LDOS at a hole-rich 
site p % « 0.30; we see that A PG {T = 40K) « 24meV 
and Ad(T = A0K) = 7.1meV (marked by arrows in the 
plots) and they both vanish at T*(i) » T c (i) « 72K. At 
a hole-poor site of p l = 0.9, A PG (T = 40K) « 37meV 
and Ad(T = A0K) = 4.5meV. Following the tempera- 
ture evolution we see that T c (i) ~ 72K that is also the 
critical temperature T c (0.21), and Apq vanishes at much 
larger value T*(i) ss 105-fT. Thus, the pseudogap phase 
of overdoped samples is composed mainly by Ap G - 

The above results lead us to many conclusions con- 
cerning the HTSC measured properties: i) The charge 
inhomogeneities and the charge segregation potential oc- 
cur due to the formation of (almost zero doping) AF re- 
gions. Consequently, its effect is more intense at low dop- 
ing compounds, leading to large superconducting ampli- 
tudes Ad and Ape in underdoped samples. On the other 
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FIG. 3: The p = 0.11 LDOS. Top panel, A PG (T) at a hole- 
poor puddle (pi = 0.021). A PG (T = 40K) « 80meV. The 
superconducting gaps Ad(T = 40K) = 12meV are marked 
by arrows. Below, the set of LDOS curves at a hole-rich 
grain, (p(i) w 0.20) with A PG (T = iOK) w 50meV and with 
A d (T = 0) = 33meV. In the inset A d (i, T)xT for each case. 
Both A PG (i) and A d (i) vanish near T = 147K. 



hand, at far overdoped region, the large Coulomb repul- 
sion in the hole-rich grains destroys the EPS transition 
and the superconductivity, it) As in granular supercon- 
ductors, the resistivity transition occurs by Josephson 
coupling among the intragrain superconducting regions. 
Hi) In this scenario, the pseudogap is due to the weakly 
localized energy levels in the two dimensional puddles 
and it is not directly related to the intragrain supercon- 
ductivity, although both phenomena are originated by 
the segregation potential V 9 &. This different process was 
demonstrated by the distinct behavior of the two signal 
with applied magnetic fields and also with the tempera- 
ture in tunneling experiments 27j • iv) The observed dif- 
ference in the LDOS form, called "coherent" and "zero 
temperature pseudogap" @ , is due to the LDOS at hole- 
rich and hole-poor locations respectively. This distinc- 
tion is clearly seen in the overdoped calculations (Fig.fllJ)) 
where most LDOS presents coherent peaks in opposi- 
tion to the underdoped case (as in Fig.©) where the 
rounded and ill-defined peaks are more abundant, v) 
The hole-poor regions have basically one electron per 
unit cell and the large Coulomb repulsion generates a 
large asymmetry between electron extraction and injec- 
tion as measured by the STM experiments 04^1 ■ The 
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FIG. 4: LDOS for overdoped p = 0.21. At the top, A PG {T) 
at a hole- rich grain (p; » 0.30) yielding A PG {T = 40K) « 
24meV and A d (T = 407\) = 7meV marked by the ar- 
rows. Below, Apc{T) at an "insulator" grain (p(i) ~ 0.09), 
A PG (40K) w 37meV and A d {T = 40.ff) = i.5meV. In the 
inset, we show also that Ad(T) x T at the two locations. 
It is important to notice that Apg(i) remains even above 
T = 100K>> T c (i) for some hole-poor grains i. 



LDOS at hole-rich regions, like a homogeneous system, 
are expected to be more symmetric, vi) The observed 
anomaly or kink measured in the LDOS at very low 
bias[a-Q was demonstrated to be due to the low tem- 
perature superconducting gap A«j. Fig.© shows clearly 
a small kink marked by the arrows. In general, it is 
more visible in underdoped samples because the Ape is 



much larger than A<j. vii) For underdoped samples, as 
in Fig.©, Ad and Ape remain finite much above and 
vary very little around T c (p) ~ 65K. The presence of 
the superconducting gap and its quasiparticle dispersion 
above T c (p) was measured in weakly underdoped Bi2212 
by ARPES experiments [12| that also showed that the 
gap (Ape) almost does not change around T c (j>). viii) 
The ARPES experiment of Lee et al[10] has measured 
increasing gapless Fermi arcs along the nodal (tt, tt) re- 
gion with increasing doping above T c {p). As we have 
shown, the d-wave A^ remains above T c (p) in the un- 
derdoped regions but tends to vanish close to T c (p) for 
overdoped compounds. Assuming that the weakly bound 
states (Apg) are due to the random phase boundary po- 
tential Vgb-, and that they occur mainly along the Cu — O 
(antinode) direction, together with the d-wave behavior 
of Ad is agreement with the increase of the gapless Fermi 
arcs with p above T c (p)[10}. ix) The calculations show 
that, in general, for any compound, the pseudogap Ape 
is smaller but with well defined peaks at the hole-rich 
regions than in the hole-poor ones. Consequently, the 
measured larger Apc{i) gap values located at the low 
density grains with local (AF) insulator behavior have 
lower local conductivity (dl/dV) than the smaller gaps 
(and higher densities) as verified by Pasupathy et alQ- 
x) The STM measured relation 2A/K B T*{i) w 8.00] is 
reproduced closely by the optimal and the p — 0.21 sam- 
ple, considering A = Apc(i) and its associated T*(i). 

In summary, we have proposed an EPS transition to 
describe the inhomogeneous charge distribution of HTSC 
generated by the lower free energy of the (undoped) AF 
domains. This approach yields potential wells with shal- 
low bound states that reduce the kinetic energy and fa- 
vor the superconducting pairing. The calculations obtain 
the inverted bell shape critical line T c (p) x p, distinguish 
clearly the LDOS pseudogap Ape from the intragrain 
superconducting gap energy Ad and provide simple phys- 
ical interpretations to many different STM and ARPES 
results. 
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